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摘 要 基于 弱 非 线性 热 声 理论 , 对 热 声 换 热 属 的 基本 换 热 特性 进行 了 理论 分 析 . 在 前 人 研究 基础 上 , 增加 了 换 热 恬 内 纵 
向 温度 梯度 的 影响 ， 获 得 了 二 阶 周 期 平均 热流 的 解析 表达 式 ， 并 定义 了 适合 表征 热 声 换 热 句 换 热 特征 的 无 量 纲 换 热量 . 计 
算 分 析 了 无 量 纲 参数 声 导 率 Ya, 对 流 与 压缩 性 比 ,动态 雷诺 数 Rew 对 换 热 器 时 均 换 热 量 的 影响 . 分 析 表 明 ,， 换 热 器 产 
生 时 均 吸 热 还 是 时 均 放 热 特性 主要 由 临界 声 导 率 幅 值 决定 , 存在 时 均 纵 向 温度 梯度 时 , 还 取决 于 Rew; 声 导 率 的 相位 角 对 
时 均 热 流 的 影响 相对 较 小 , 但 其 影响 会 随 着 |T| 和 |Yal| 的 增 大 而 增 大 ; |T| 的 存在 增强 了 换 热 ， 存在 临界 的 Rew 值 使 得 
二 阶 周 期 平均 热流 取得 极 值 . 本 文 所 定义 的 无 量 纲 换 热量 能 合理 表征 换 热 器 换 热 特性 ,对 热 声 换 热 避 基本 换 热 特征 的 研究 
将 为 进一步 考察 其 实际 复杂 工作 机 理 提供 理论 指导 . 
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Abstract Based on the veal oh 
heat transfer characteristics Se heat exchanger are analyzed in this paper. Mean time 
longitudinal i 


near theory of thermoacoustic energy conversion, the basic 


bi is considered. A nondimensional heat flux is defined and it’s mainly 
influenced by Pet S 和 as dynamic Reynolds number Rew，acoustic admittance Ya, and the 
temperature ratio 三 . 二 analytical expressions Show that there exists critical acoustic admittance 


making the heat flux being zero. And the existence of |T'| enhamced heat transfer. This work can 


improve the understanding of heat transfer in oscillating flow 二 can be useful for intensive study of 
0 practical thermoacoustic heat exchangers. XK ” 
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0 引 言 «x 


热 声 热机 是 基于 热 声效 应 工作 的 ， 因 其 无 运动 
部 件 、 结 构 简 单 、 工 质 无 污染 、 可 靠 性 高 、 能 利用 低 
品位 能 源 等 优点 而 受到 了 很 多 研究 者 的 重视 . 热 声 
热机 的 关键 热力 部 件 包括 换 热 器 和 回 热 器 ， 其 中 换 
热 吉 使 回 热 需 的 两 端 形 成 温度 梯度 ， 实 现 系 统 与 外 
界 的 热量 交换 ,以 维持 热 声效 应 . 然而 , 由 于 交 变 流 
动 本 身 的 复杂 性 和 多 样 性 ， 至 今 我 们 对 交 变 流动 与 
换 热 规律 的 认识 与 稳 态 换 热 相 比 仍 相差 其 远 。 

热 声 换 热 器 由 两 部 分 组 成 ， 一 部 分 是 外 界 热 源 
用 以 提供 换 热 器 所 需 的 外 部 工作 环境 ， 一 部 分 是 系 
统 内 部 的 交 变 流动 工作 流体 ， 换 热 器 的 换 热 效果 由 
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这 两 部 分 看 合 决定 。 本 文 将 针对 内 部 流体 的 换 热 进 
行 理论 分 析 . 内 部 工作 流体 多 采用 惰性 气体 ， 在 热 
声 系 统 工作 过 程 中 ， 所 有 状态 参数 与 流动 参数 都 是 
波动 的 ， 换 热 过 程 尤为 复杂 -3 热 声 换 热 器 中 的 
热量 交换 是 由 动态 换 热 和 二 阶 的 周期 平均 换 热 两 部 
do eee Gt 

， 不 能 直接 表征 换 热 器 的 换 热 性 能 . 在 热 声 系 
我 们 更 关心 的 是 周期 平均 能 量 转 换 效 应 , 对 二 
阶 的 时 均 换 热 研究 属于 非 线性 热 声 动力 学 范畴 . 通 
过 对 热 声 系统 的 能 流 分 析 可 知 回 : 不 论 声波 振幅 大 
小 , 热能 与 声 能 转换 本 质 上 是 非 线性 的 , 而 热 声 转换 
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正 是 依 束 于 这 种 非 线性 效应 而 工作 的 ， 即 对 于 热 声 
系统 而 言 ,能 量 转换 必须 要 考虑 二 阶 周 期 平均 效应 ， 
基于 此 ， 已 有 研究 中 67 采用 的 理论 分 析 模 型 将 换 
热 器 视 为 无 限 长 等 温 流 道 ， 并 着 重 分 析 了 声 导 率 对 
换 热 器 的 影响 但是， 实际 换 热 器 的 长 度 只 有 约 两 
倍 当地 流体 振东 位移 幅 度 ， 存 在 显著 的 时 均 温度 分 
布 以 及 进出 口 效应 .为 简单 起 见 ， 本 文 首先 考察 时 
均 温度 梯度 的 影响 ， 对 换 热 器 基本 换 热 特性 做 更 深 
一 步 的 考察 ， 进 出 口 效应 等 强 非 线性 问题 的 研究 将 
作为 后 续 进 一 步 的 研究 内 容 . 

本 文 以 弱 非 线性 热 声 理论 为 基础 ， 从 原始 的 能 
量 守 便 方 程 出 发 ， 获 得 了 考虑 纵向 温度 梯度 影响 
的 平行 平板 换 热 器 的 二 阶 周期 平均 热流 的 解析 表达 
式 ， 并 对 其 进行 了 无 量 纲 化 ， 主 要 考察 了 动态 雷诺 
数 Rev， 当 地 声 导 率 Ya 以 及 对 流 与 压缩 性 比 
对 无 量 纲 二 阶 周期 时 均 热流 Go 的 影响 ， 分 析 获得 
了 热 声 换 热 器 的 一 些 基 本 换 热 特性 . 与 以 往 采 用 稳 
态 分 析 方 法 引入 努 塞 尔 数 定义 来 强 0 
与 时 均 温 差 不 同 ， 本 分 析 方法 基于 基本 的 热 声 i 
0 
符 性 


1 二 阶 周期 平均 热流 的 sh 


对 能 基 守 恒 方 程 的 微分 焉 了 多 分 与 时 
均 运 算 ， 获得 含有 壁面 换 热量 的 时 均 能 量 方程 如 下 


所 示 : 


(te) [to =0 


其 中 ，z 为 声 传播 方向 ，(pVih)o 代表 z i 
ns (et) 代表 z 方向 的 动能 流 ， 为 三 阶 量 ， 


2 
0 

可 忽略 (常用 换 热 器 中 的 速度 幅 值 一 般 在 10 m/s 量 

级 , 相 较 于 时 均 灼 流 可 忽略 , 如 速度 达到 100 m/s 量 

级 , 或 变 苓 面 附近 则 需要 谨慎 处 理 ); (过 .六 ) 指 


共性 力 做 功 ， 为 三 阶 基 ,可 忽略 ;( -5 ) 代表 
方向 的 导热 ,气体 导热 系数 冬 小 , 与 时 均 炊 流 的 变化 


相 比 一 般 可 忽略 ; Qo 代表 流体 壁面 的 时 均 热 流 . 因 
此 ， 对 于 换 热 器 ， 流 体 壁 面 的 时 均 热 流 就 近似 等 于 
流 道 进出 口 的 时 均 灼 流 差 ， 

根据 线性 热 声 理论 ,无限 长 流 道内 存在 时 均 温 
度 梯度 时 的 速度 波动 、 压 力 波动 以 及 温度 波动 都 可 
获得 解析 解 四, 将 这 些 解 析 表 达 式 代入 式 (1) 中 , 妈 
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流 周期 平均 值 保留 到 二 阶 精 度 ， 忽 略 高 阶 项 ， 得 到 


[4,8,9] ， 


i 局 fi— fe 
Go= neB {eT (+ 中 


poCplUnl? 
3A01 — or) 1 — fp 


其 中 ，P 表示 压力 ,UV 表示 体积 流 率 ， 上 标 * 表 
示 共 罗 ， 下 标 “1” 表 示 一 阶 波动 量 ，c” 为 普 朗 特 
数 ，Re 表示 取 实 部 ，Im 表示 取 虚 部 运算 。 对 于 
平板 结构 ， 截 面 平均 的 黏 性 函数 f, = tanh(2)/2, 
2 = (1 十 人 VRew/32; 截面 平均 的 热 函 数 f= 
2Z). 热 、 黏 函数 表征 的 是 参数 的 横向 分 
是 截面 形状 以 及 水 力 直 径 与 黏 性 穿 透 深 
态 雷 ed 0 0 


mf tof) Pe)} © 


tanh(oZ)/(o 
布 的 影 响 ， 


得 热 、 ee oe ed 
梯度 , 但 这 些 效应 相对 对 流 项 的 影响 较 小 , 因此 , 为 
了 简便 起 见 , 暂时 忽略 换 热 器 中 的 温度 分 布 对 热 、 生 
函数 及 其 温度 非 线性 分 布 的 影响 ， 只 考虑 压力 波动 
与 体积 流 率 波动 的 分 布 特性 . 将 压力 梯度 与 体积 流 
率 梯度 均 转 换 为 压力 与 体积 流 率 的 表达 式 ， 从 而 获 
得 时 均 热流 的 表达 式 如 下 : 


» Aw|P 办 后 
NI 


eb 


DF 一 六 /poaoUiN 

加 人 Ti (2 
站 

;\ i |) + 


eixaro-Di (Se))+ 


poaoU1 Pe Cp dion, 一 — 1)(f: — = | 


AP aow dz (1 — f,)(1—o0) 
全 Oo = 452 oo 得 到 无 量 网 二 阶 周期 平均 
热流 ， 
i fr— fr 
HS Re 人 (1 )* 
, je 
fit 0 -D+ Ye 


i ll 2Thn(fe+of’) 
re) [bP os)* 


[re x (1+(y—1fr)Ya)+ 


其 中 , 无 量 纲 参数 Ya = oo 了 = CS 
AP aow dx 

由 式 (4) 可 以 看 出 6@o 是 流体 的 物性 参数 7、Pr、 
流动 参数 Rew、 声 导 率 Ya 以 及 对 流 与 压缩 性 比 卫 
的 函数 . 卫 = 0 即 不 考虑 温度 梯度 影响 时 , 式 (4) 与 
文献 57] 中 时 均 换 热 表 达 式 有 相同 形式 . 热 声 热机 
中 的 常用 工 质 为 氨 气 与 氮气 ,本文 以 氧气 为 例 ， 取 
Pr = 0.67、 7 二 1.67, 着 重 分 析 Rew、Ya 及 本 对 90 
的 影响 , 即 6o = f(Ya, T, Rew). 


2 二 阶 周 期 平均 热流 的 影响 因素 分 析 

在 热 声 系 统 中 ， 换 热 器 的 横向 尺度 b (对 于 平板 
结构 为 流 道 半 宽 ) 应 与 热 穿 透 深度 5 相当 或 更 小 , 以 
保证 工作 气体 的 热 交 换 . 但 由 于 换 热 器 通道 在 实际 加 


工 过 程 中 存在 的 困难 同时 要 考虑 过 窄 流 道 带 来 的 限 代 和。 -2 和 有 ee 


增加 问题 , 6/64 一 般 被 控制 在 1~1.5 之 间 . 由 
看 出 ,常温 条 件 下 氮气 的 Rew 取信 范围 是 482107. 


小 


表 1 ”氢气 的 动态 雷诺 数 变化 放 
Table 1 The value of dynami nolds number 


as a functi Ok 
b/6k 1.0 1.1 1.2 1.3 1.4 1.8 
Rew 48 58 69 81 94 107 


利用 DeltaEC 软件 分 别 模拟 计算 了 驻 波 发 动机 
和 行 波 发 动机 两 种 常用 的 热 声 热机 系统 ， 得 到 了 换 
热 器 处 的 声 导 率 ,结合 已 有 资料 中， 计算 时 |Ya| 的 


取 值 范围 扩大 到 : 0.001~4. |Ya| = 1 意味 着 局 部 声 。 
阻抗 幅 值 等 于 poco， 即 大 空 管 行 波 特 征 阻抗 值 _XAO 


在 实际 加 热 器 中 ， 由 于 回 热 器 与 热 缓冲 恩 中 
冷 流体 影响 ， 换 热 器 纵向 温度 分 布 会 不 均匀 ,温差 
可 达 几 十 K, 其 典型 尺寸 为 50 mm, 卫 能 达到 10 的 
量 级 , 因此 分 析 中 我 们 取 研 的 范围 为 : -20 ~ 20. 
2.1 动态 雷诺 数 对 无 量 纲 二 阶 周期 平均 热流 的 影响 

图 1 给 出 了 无 量 纲 二 阶 周 期 平均 热流 随 动态 雷 
庄 数 的 变化 规律 ， 温 度 梯度 为 零 时 如 图 1(a) 所 示 ， 
除 临界 |Yal 值 附近 , 不 同 的 |Ya| 值 下 , Bo 随 Rew 
具有 相似 的 变化 规律 ， 即 Rew 一 0 与 Rew 一 co 
时 ，@o 分 别 趋 于 0 与 一 个 较 小 的 有 限 值 ， 并 在 有 限 
Rew 时 存在 一 个 峰值 ; 不 同 的 |Ya| 值 有 着 不 同 的 峰 
值 60 及 对 应 的 Rew, 该 Rew 值 在 常规 换 热 需 的 取 
值 范围 内 ; 同时 9o 的 正 负 号 即 吸 热 或 放 热 也 取决 于 
IYal 值 的 大 小 , 即 存在 临界 |Ya| 值 ( 约 2.1 左右 ) 使 
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得 换 热 量 为 零 . 

图 1(b) 为 存在 温度 梯度 的 情况 ， 极 小 与 极 大 
Rew 时 ，@o 仍然 趋 于 零 和 较 小 有 限 值 ; 但 此 时 6@o 
随 Rew 的 变化 会 出 现 两 个 反 向 峰值 , 不 存在 |Ya| 值 
使 得 6o 沿 Rew 完全 处 于 吸 热 或 者 放 热 状态 ， 存 在 
临界 的 Rrw 值 使 得 @o 经 过 零点 ; |Ya| 值 越 大 对 应 
的 临界 Rev 越 大 . 实际 换 热 器 中 可 以 通过 优化 流体 
的 流动 参数 达到 最 佳 的 换 热 效果 . 


< 


1000 


Re, 
(a) T=0, $B(Ya)=60° 


0.01 0.1 1 10 100 
Re, 
(b) [=2, B(Ya)=60° 


图 1 无 量 纲 二 阶 周 期 平均 热流 随 动态 雷诺 数 的 变化 规律 
Fig. 1 dimensionless second-order time-averaged heat fux vs 


Reynolds number 


2.2 声 导 率 对 二 阶 周期 平均 热流 的 影响 规律 

以 下 分 别 考察 声 导 率 的 幅 值 和 相位 角 对 二 阶 周 
期 平均 热流 的 影响 规律 . 图 2 给 出 了 无 量 纲 二 阶 周 
期 平均 热流 随 声 导 率 的 幅 值 的 变化 规律 ， 

图 2(a) 表明 , 存在 临界 的 |Ya| 使 得 Oo 经 过 零 
点 , 不同 的 Rew 值 对 应 的 临界 值 不 同 ， Rew 值 越 大 
对 应 的 临界 | 如 | 值 越 大 , 但 该 数值 变化 不 大 ， 所 考 
察 Rew 值 范围 内 临界 |Ya| 值 均 在 2.1 左右 ; 图 2(b) 
表明 , 考虑 温度 梯度 时 ，6@o 随 |Ya| 的 变化 存在 一 个 
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先 减 小 后 增 大 的 过 程 ， 且 存在 临界 的 声 导 率 的 幅 值 
使 得 二 阶 周期 平均 热流 为 零 ， 声 导 率 的 幅 值 小 于 临 
界 值 时 ,二 阶 周 期 平均 热 沈 大 于 零 , 此 时 为 放 热 句 ， 
反之 为 吸 热 带 ; 记 值 不 等 于 零 时 , 随 着 |Ya| 增 大 , 换 
热 句 成 为 放 热 器 ， 且 换 热 效果 显著 增加 . 


| 玛 | 
(a) T=0, B(Ya)=60° 


|Yal 
(b) T=5, B(Ya)=60° 


图 2 无 量 纲 二 阶 周期 平均 热流 随 声 导 率 的 幅 值 的 变化 规律 


Fig. 2 dimensionless second-order time-averaged heat flux 9 


«~ 
图 3 给 出 了 无 基 纲 二 阶 周期 平均 热流 随 声 导 率 
的 相 角 的 变化 规律 , 存在 温度 梯度 时 相位 角 才 对 Go 
有 影响 ，@o 随 5(Ya) 正 弱 变 化 。 随 着 有 和 |Yal 
的 增 大 ， 二 阶 周 期 平均 热流 对 声 导 率 的 相位 角 敏 感 
度 均 增加 ， 即 相位 角 对 二 阶 周期 平均 热流 的 影响 增 
大 . 6o 随 B(Ya) 的 变化 在 6(Ya) = 一 10° 时 存在 平 
衡 点 , 在 86(Ya) = 一 100° 或 8(Ya) = 80° 时 分 别 取 
得 最 大 或 最 小 值 . 在 实际 换 热 器 中 ， 通 过 调整 声场 
条 件 可 以 获得 更 合适 的 换 热 效果 . 
2.3 无 量 纲 参数 三 对 二 阶 周期 平均 热流 的 影响 规律 
图 4 给 出 了 无 量 纲 二 阶 周期 平均 热流 9o 随 无 
基 纲 参数 荆 的 变化 规律 , 正人 负 温 度 梯度 下 ,eo 的 变 
化 具有 对 称 性 . 图 4(a) 表明 , 不 同 的 |Ya| 值 下 eo 


acoustical admittance 
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随 荆 的 变化 规律 相似 ，eo 随 着 |T| 值 的 增 大 而 增 
大 , 即 换 热 效果 越 来 越 好 , | Ya| 值 决定 了 9o 随地 变 
化 的 敏感 度 ; 图 4(b) 表明 , 不 同 的 Rew 值 下 eo 随 
的 变化 规律 相似 ，9o 随 着 |T| 值 的 增 大 而 增 大 ， 
即 换 热 效果 越 来 越 好 . 
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图 3 无 量 纲 二 阶 周期 平均 热流 随 声 导 率 的 相 角 的 变化 规律 
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(b) |Yal=1, B(Ya)=60° 
图 4 无 量 纲 二 阶 周期 平均 热流 随 无 量 纲 参数 忆 的 变化 规 得 


Fig. 4 dimensionless second-order time-averaged 
heat flux vs T 


ay 


3 结 论 


本 文 以 弱 非 线性 热 声 理论 为 基础 ， 研 究 在 小 振 
幅 波动 条 件 下 换 热 器 的 二 阶 时 均 换 热 特 性 ， 获 得 了 


| 了 平板 通道 内 二 阶 疝 引 
平均 热流 的 解析 式 ， 给 出 了 交 变 流动 换 热 属 的 热 冰 


并 指 出 存在 由 Rew 值 俩 得 


得 二 阶 膨 期 平均 热流 为 堆 ， 当 实际 
于 临界 值 时 ， 振 获 流 体 从 外 热源 能 部) 为 吸 热 器 ; 当 
实际 声 导 率 的 幅 值 小 于 临界 低 暑 、 振荡 流 体 向 外 放 
热 , 为 放 热 器 ; 存在 温度 梯度 W], 随 着 |T| 和 |Ya| 的 
增 大 , 声 导 率 的 相 角 对 Gu 的 影响 越 大 ; 另外 考察 了 
无 量 纲 参数 对 二 阶 周期 平均 热流 的 影响 ,结果 表 
明 ，6u 随 着 |T| 值 的 增 大 而 增 大 . 

以 上 研究 为 进一步 考察 交 变 流动 换 热 器 的 换 热 
特性 提供 了 依据 ， 对 提升 换 热 器 的 工作 性 能 及 改进 
换 热 器 的 结构 设计 具有 重要 的 指导 意义 x 
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